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Sum rules are used to calculate the corrections to the hydrodynamic collective frequencies of a 
harmonically trapped Bose gas at T — 0, due to finite size efi'ects. We show that, with logarithmic 
accuracy, the relative corrections behave like {aho/R)'^^og{R/aho), where R is the Thomas-Fermi 
radius of the condensate and aho is the oscillator length. Results are given for both spherical and 
axi-symmetric trapping. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj, 67.90.+Z 



I. INTRODUCTION 



'■^ I After the experimental realization of Bosc-Einstein condensation (BEC) in atomic trapped alkali gases ||l[ the study 

fT^ of their dynamical properties has become an important subject of experimental and theoretical research (see M and 
O^ ' references therein). The measurements of the collective frequencies have confirmed the validity of the mean field 
approach l2l-H] in describing the dynamical behaviour of a trapped Bose-Einstein condensed gas at zero temperature, 
properly accounting for the important role played by tvifo-body interaction and by Bose statistics. The high precision 
of these frequency measurements is stimulating new perspectives in the study of finer effects like the temperature 
dependence of the collective frequencies, the occurrence of beyond mean field corrections [g|, etc. In this perspective 
^h ] it is important to have a full analytical control of the values of the collective frequencies predicted by the mean field 
^^ ■ scheme at T = 0, including finite size corrections which are neglected in the Thomas-Fermi limit. The purpose of the 
^D ■ present work is to calculate these corrections in the framework of the mean field Gross-Pitaevskii theory, employing a 
J^ ,; sum rule approach p]. In general sum rules provide only an upper bound to the frequency of the lowest modes. In this 
^^ ' paper we will show that, due to the harmonic nature of the confinement, the first corrections to the Thomas- Fermi 
limit can be calculated with logarithmic accuracy using this method. We will show that the corrections are fixed by 
the kinetic energy term which, in these trapped Bose-Einstein condensed systems, plays the role of an effective surface 
tension 0. 

O . II. SUM RULES AND COLLECTIVE MODES 

O 

The usefulness of the sum rule technique in studying the dynamical behaviour of trapped Bose-Einstein condensates 
k> _ has been already discussed in several papers (see, for example, PH)- This method, typical of many-body physics, is 
'^j ■ based on the calculation of the energy-weighted moments rup — fdu) S3^{uj) {fujj)^ of the dynamic structure function 
C^ ' associated with a given operator J^ . The dynamic structure function, at zero temperature, takes the form 

5^(^) = ^|(n|.F|0)p5(fic^-fi^„o), (1) 

n 

where fiw„o — fi^n — Swo is the excitation energy of the eigenstate \n) of the Hamiltonian, whose ground state is |0). 
At T = the quantity Sj^(uj) is related to the dynamic response function of the system relative to JF through the 
expression H 
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In this work we make the following choice for the Hamiltonian 
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where the parameter g = Anfi^a/M is the strength of the two-body contact interaction, fixed by the s-wave scattering 
length a. Ratios between different sum rules provide rigorous upper bounds 



huj 
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(4) 



to the energy of the lowest state excited by the operator J-'. These bounds coincide with wiowcst only if the strength 
of the operator J^, defined by Eq. (||), is exhausted by a single mode. 

By using the closure relationship involving the eigenstates |n), the energy- weighted and the cubic moments can be 
easily rewritten in the form of sum rules as follows 



mi 



ma 
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(5) 
(6) 



where the averages are taken over the ground state of the system which, in the mean field scheme, is provided by the 
stationary solution of the Gross-Pitaevskii equation ||l^,Q . 



Let us first consider the simplest situation where the harmonic confinement is isotropic, i.e. ui^ 



^h< 



In this case the calculation of the above moments is straightforward (see Appendix ^ if the operator T coincides 
with the monopole 



TV 
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where rf 



zf, or quadrupole 
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operators, yielding the following results for the ratio \frn%fm\ 
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-Bkin = (X]i=iPf)/2^^ and Eho — -^'-^ho(Si=i ^?)/2 are the expectation values of the kinetic and 
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InEqs. (!)-( 

harmonic oscillator energies respectively, -EkinJ_ = {J2i'=i(Pxi +Pli))/'^^ = 2i?kin/3 and £^ho_L — M Ljl^{J2i=i{xl 
))/2 = 2 -Eho/3 are their corresponding radial contributions. In deriving Eq. M) we have also used the virial identity 

2£;ki„ - 2Eho + 3Eint - (11) 

relating in an exact way the various contributions to the ground state energy 

E — -Ekin + -E'ho + Eint (12) 

where Eint = 5(X]i<7" ^('^^ ~ '^j)) ^^ ^^^ average value of the mean field interaction energy. 

Another possibility for calculating the monopole and the quadrupole frequencies makes use of the ratio y^mi/mZi , 
which involves the inverse energy weighted moment. In general the following inequality holds 



\fm^JmZ[ < y^ms/mi (13) 

showing, by comparison with (Q), that the lowest excitation frequency is in general better approximated by the 
ratio ■\/mi/m,_i than by y^m^/mi. Differently from nii and 7713, the moment m_i cannot be written in terms of 



commutators, but is related to the static response function via the following general relationship (compressibility sum 
rule §) 



X^(0) = -2^ 
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(14) 



which provides a natural rule for calculating ra-i through the explicit determination of the polarization {T + ^^) 
^X^ induced in the system by an external static field interacting through the Hamiltonian ^{J- + J-^). 
We present in Appendix H the details of the derivation and report here only the final results 0Jl3l : 



^Mi,-i — 2uJho 
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(15) 
(16) 



are, respectively, the deformation of the confining 
|) and dig ) show that the monopole frequency ujmi _i 

N 



where s = {ul — LUy)/{ujl + cOy) and S — (y^ 



potential and of the condensate in the x-y plane. Equations 

is determined by the A^-dependence of the square radius (r^) — (X]i=i ^1)/^^ while the quadrupole by the dependence 
of the deformation 5 of the cloud on the deformation e of the trap. In (llfl) we have explicitly taken the limit of a 
symmetric trap. 

It is worth pointing out that the results (]9|)-(pX|) and (|l^)-(16) are valid in general in the framework of the mean 
field Gross-Pitaevskii theory and apply, in particular, to both repulsive and attractive forces Q. In particular 
they reproduce the correct values for the frequencies in both the relevant limits of the non interacting gas, where 
-Ekin = ^ho, and of the Thomas- Fermi approximation where the kinetic energy is negligible with respect to i?ho pq |. In 
the former case both the monopole and quadrupole frequencies approach the value 2a;ho, while in the latter they tend 
to the analytic hydrodynamic values v5t^ho and v^i-^ho for the monopole and quadrupole oscillations respectively, as 
predicted in Ref . |3| . This feature allows us to conclude that in the hydrodynamic regime a single mode exhausts the 
sum rules 7743, mi and to_i. In the following section we will show that, with logarithmic accuracy, the same is true 
also if one includes the leading finite size corrections to the Thomas-Fermi limit. 



III. FINITE-SIZE CORRECTIONS IN SPHERICALLY SYMMETRIC SYSTEMS 

In the Thomas- Fermi regime the asymptotic behaviour of the kinetic energy per particle of a harmonically isotropic 
Bose gas at T = is given by [[l6|jl7|] 
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where 
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(18) 



is the Thomas-Fermi radius of the condensate, with aho = \/^JMuj\^, and C ~ 1.3 is a numerical factor. Such a 
behaviour arises from a careful analysis of the shape of the order parameter near the border of the condensate, where 
the Thomas-Fermi approximation fails and the wave function of the condensate should be calculated with proper 
accuracy. By using the virial identity (|l^), and the usual definition /i = dE/dN for the chemical potential one can 
easily show that the first corrections to E\^q/N and Eint/N beyond the Thomas-Fermi approximation are given by: 

(19) 

(20) 
where 
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is the Thomas-Fermi chemical potentiaL Note that in deriving the above equations we have consistently taken into 
account all the terms of order 1/i?^, which renormalize the constant C inside the log argument. Eqs. (|l^)-(|2y) were 



also derived in Rcf. fl^^ employing a boundary layer theory. 

t the first corrections i 

-Ekin 35 /aho\'* , / R 



a. Rcf. 113 by ^ 
From Eqs. (|l7| ) and ( [19|) one finds that the first corrections to the ratio i?kin/£'ho can be written as 



m'--c^ '-' 



Eho 3 V i? / y CcL 

and becomes smaller and smaller in the large i?, Thomas-Fermi limit. Equations (M) and (ll^) also reveal that the 
kinetic energy plays the role of an effective surface energy, providing the first correction to the total energy (|l2| ) in 
the Thomas- Fermi regime. 

By using the above results it is immediate to determine the finite-size corrections to the hydrodynamic value of the 
frequencies predicted by Eqs. (O) and (UQ). It is interesting also to evaluate these corrections starting from Eqs. dlq) 
and (|l^). Appendix y is devoted to the technical details of this calculation, which shows that ratios uj^^i and uji-i 
do coincide for both the monopole and the quadrupole modes, apart from a different value of the coefficient inside 
the log argument. The results for the collective frequencies, including the first corrections due to finite size effects, 
take the form 
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lim UOQ — V 2 Who 
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(23) 

(24) 



with Cm — C e^^l^ and Cq = C e^l^ for cji._i, while Cm — Cq = C if one uses the estimate loj,^\. Notice that 
our results differ from the ones found in Rcf. |fj], where corrections behaving like (oho/^)"' (without the logarithmic 
factor) were obtained by employing a different procedure based on a perturbative solution for the excited state wave 
function. 

The physical role played by the kinetic energy and the R dependence of the finite size corrections can be understood 
in terms of the classical picture developed in Ref. |0] for the surface modes of a Bose condensed gas confined by a linear 
potential. In this case the resulting dispersion relation takes the form uP' — Fq/m plus q^ log q corrections originating 
from the kinetic energy term, where q is the wave vector of the surface wave. In this equation F — mujf^^R is the 
harmonic force evaluated at the Thomas-Fermi radius R. The correspondence with our result ( |24| ) for the surface 
quadrupole mode is achieved by imposing the discreteness condition q ~ \/ R to the wave vector. 

It is also worth noticing the different sign of the finite size corrections entering the monopole and quadrupole 
frequencies. Furthermore the correction for the quadrupole is about 5 times larger than for the monopole. 

It is finally interesting to compare the finite size corrections predicted by Eqs. (|23| ) and (gj) with the ones due to 
beyond mean field effects [Q. These are fixed by the gas parameter 

1^2/5 / „ \12/5 

aM^) = ^[N"' — ) (25) 

»7r \ Oho / 

calculated at the central density n(0) of the gas. These effects are not included in the mean field Gross-Pitaevskii 
theory and take into account the changes in the equation of state originating from quantum correlation effects according 
to the theory of Lee- Yang-Huang |]l8|. The dependence of the gas parameter (^)on the relevant physical variables 
(iV, a and Oho) is quite different with respect to the one exhibited by the ratio (22|). Furthermore in Ref. [g[ it was 
found that only the monopole oscillation is affected by beyond mean field effects. This should introduce enough 
flexibility to distinguish, in the study of the collective frequencies, between finite size and beyond mean field effects. 
First experimental efforts in this direction were made on *^Rb where the value of the scattering length can be tuned, 
thereby amplifying the effects of quantum correlations [Q . 

IV. RESULTS FOR DEFORMED TRAPS 

Until now our analysis was restricted to the case of spherical traps. It is interesting, especially from the experimental 
viewpoint, to extend the above results to the case of deformed external potentials. We will focus on harmonic traps 
having axial symmetry [ujx ^ LUy = lu±_ and ujz = Acjj^, the parameter A characterizing the shape of the trap). Having 
proven in the previous section that the accuracy of the sum rule approach is logarithmic, we will present in the 
following only the corrections to the estimate y^m^/rrii, whose calculation is straightforward. 



We are interested in particular in the study of the collective excitations described by the 
operator (g) and by the operator 



N 



^-E 



-pzf 



2, m — 2) quadrupole 



(26) 



which accounts for the coupling between the i — 2, m — and the the i = excitations. The parameters a and /? fix 
the relative weight of the two modes. Moreover one can show that by treating them as variational parameters within 
a sum rule calculation, one recovers the correct dispersion for the two decoupled modes in Thomas-Fermi regime, first 
obtained in Ref. ||^ by solving the linearized hydrodynamic equations of superfluids. 

Let us first discuss the m = 2 mode. By using the results (A3) and (A4) reported in Appendix^ it is straightforward 
to find the result 
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where the Thomas-Fermi chemical potential is the natural generalization of Eq. (Ell) 



AiTF = ^MujIrI , 



with the radial size of the condensate given by 
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(28) 



(29) 



and a± = ^Jh/Mu)^. 

For the m = mode, by using Eqs. ([A5|) and (RS), we find 



u;2(m = 0) = 



1 a2(io E^, - 2 Skin) + 2a/3(10Sho^ - 2 Skin.) + /3'(6 E,,,, + 2 E^nz) 
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2a/3(z2)+/32(z2) 



(30) 



for the ratio (ftw)2 = 7713 /?tii, where Skin z = (X]i=i-Pz )/2-^ and Shoz = Ma;2(^^^j^ Zi2)/2 are the axial contributions 
to Skin and Sho respectively. By imposing that the ratio ( pO| ) be stationary with respect to the variation of a (or /?) 
one can determine the decoupled eigenfrequencies. We present here only the explicit results in the simpler cases of 
strongly deformed traps. In the limit A <C 1 (cigar) one finds that the lowest of the two frequencies of Eq. (pO) reads 



lim id (ra — Q] 
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(31) 



Conversely when A :^ 1 (disk) one gets 



lim Lo (m = 0) = 
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According to Eqs. (P7|), (31) and (p^), in order to make quantitative estimates of the finite size corrections one 
needs an explicit expression for Skin ± and Skin 2 ■ We have calculated these quantities by employing the method 
developed in Ref. [pT[. In Appendix O we report the results, which explicitly depend on the deformation of the trap. 



V. CONCLUSIONS 



In this paper we have used a sum rule approach to calculate the frequency of the monopolc and quadrupole modes 
of dilute Bose gas confined by a harmonic trap at T = 0. We have shown that with logarithmic accuracy the sum 
rules are exhausted by a single mode not only in Thomas-Fermi regime but also if one includes the first corrections 
due to the finite size of the could. We have provided an explicit expression for the frequency shift of the monopole 
and quadrupole modes due to finite size effects both for spherical and deformed traps. 
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APPENDIX A: RESULTS FOR THE ENERGY WEIGHTED AND THE CUBIC SUM RULES 

The mi and m^ sum rules can be calculated by carrying out the commutators of Eqs. (p|)-(p) with the mean field 
Hamiltonian (|3|). The results for ?7i_i are instead derived in Appendix y[ 
For the monopole operator (M) one finds: 

mf = ^-^N{r') (Al) 



-M 



M 



„„, - ^2-(4£;ki„ + 4£;ho + 9£;int) . (A2) 

For the quadrupole operator (^) the results read: 

mf = ^N{x'+y') (A3) 
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(£;ki„±+^ho±) . (A4) 



For the operator ( pq ) we find: 



mf = ^N{a^r^) + 2aP{z') + P^z^)) (A5) 

2h? 

mf = _(a2(ioSj^^ _ 2^kin) + 2a/3(f0£;ho. - 2£;ki„.) + /3'(6Sho. + 2Skin.)) ■ (A6) 

APPENDIX B: THE INVERSE ENERGY WEIGHTED SUM RULE 

The aim of this appendix is to determine the inverse energy weighted sum rule relative to the monopole and 
quadrupole operators (R) and (0) when the Hamiltonian of the system is given by Eq. (H), with ujx = (^y = oj^ = oj\^o- 
The procedure is based on the determination of the static response function through Eq. (pH). 

Monopole: in this case the problem is equivalent to finding the static response when the system is perturbed by an 
external field proportional to the monopole operator (|7|), which can be described by the Hamiltonian: 

N 



^pcrt=7E^^ (Bf) 



i=l 



where 7 is the strength of the perturbation. Notice that adding such a perturbation to the Hamiltonian (@) has simply 
the effect of renormalizing the frequency of the external confinement, which now becomes 



27 
M 



^7 = \/< + 77 ■ (B2) 



As a consequence, the average square radius will now depend on 7. 
By definition, the static polarizability is given by 



d{ 



2\ 



Xm{uj = 0) = N[^^^] . (B3) 



In order to evaluate (B2), one should exploit the 7 dependence of {r )-y in an explicit way. To this purpose it is 



convenient to write the stationary Gross-Pitaevskii equation in a suitable form, by rescaling the lengths, the energies 



and the order parameter as suggested in Ref. 



r ~ a-yTi 
E — huj-y El 

N 



(B4) 
(B5) 

(B6) 



with V^i(ri) norniahzed to 1. According to such notation the time-independent GP equation for ipiiji) reads 

[-V? +rl + ui|V'(ri)p]V'i(ri) = 2fiiMri) , (B7) 



where ui — Sir Na/a^, a^ = ^hjMuj^ and /i — hu-y^i. Consequently one can write the square radius of the system 
as 



(r^)^ = a'^Arf)^ = "7 / ^i"i ^1 l'/'i(i"i 



(B8) 



Note that {rf)-y exhibits a dependence on 7, since the Gross-Pitaevskii equation, in the form (B7), depends on 7 
through the mean field term proportional to ui. The derivative of (r^)^ with respect to 7 can then be expressed in 
terms of the deriv ativ e with respect to ui or N. 
By using Eqs. (B8) and (p3) one finally finds 
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(B9) 



which, combined with Eq. (^^, gives the relevant result 
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Quadrupole: in this case the idea is that in the presence of a harmonic potential deformed in the x-y plane (tOx "/= '-^y) 
the quadrupole operator (g) can be written also in the form 



Q 



1 [H,L,] 



where 
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(Bll) 



(B12) 



is the third component of the angular momentum operator. By using the relationship (nj) combined with the 
expression ( Bll ) for the quadrupole operator, one easily shows that m_^ can be written as 



o _1 {[L.,Q]) 

-1 2^hM{u^l-u^l)■ 



From Eq. (B13), taking the spherical limit, one immediately finds the result 
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,S _ 
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where the deformation of the could S and the one of the trap e have been introduced in Sec. Ill 



APPENDIX C: COMPARISON BETWEEN cja.i AND tui,_i 



We report here the fuU derivation of the finite size corrections relative to both nionopole and the quadrupole modes 
in the isotropic case. 

Monopole: It is immediate to derive the corrections to uj^^i by using the Eq. (^ combined with the Eq. ([2^). One 
finds the result 



UJm 



3 1 = VdUJho 
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It is quite simple also to calculate the finite size corrections to wi.-i starting from Eq. (|l5|). One has to use the 
identity 



Eho 



Mu^L N ' 



(C2) 



with the harmonic oscillator energy given by Eq. (|l9|). Note that the A^-dependence of the mean square radius (r^) 
originates from both the Thomas- Fermi chemical potential p^) and from the correction behaving like 1/R^ log(i?/aho), 
where R is given in Eq. (18). By carrying out carefully the calculation and taking into account all the terms of order 
1/R'^ , one finally finds 
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which coincides with the result (CI) for lu^^i within logarithmic accuracy. Note that the inequality ( |13[ ) is correctly 
fulfilled by Eqs. (CI) and (C3), the relative difference between the two estimates becoming smaller than 10^** for 

R>6a\io- 

Quadrupole: the calculation of the corrections starting from the ratio y^ms/mi is immediate if one combines Eqs. 
(|l^) and dl). This yields the result 
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The calculation of the corrections to wi^_i is less trivial since it requires the explicit knowledge of the dependence 
of the condensate deformation S on e. To this purpose we determine the finite size corrections to the square radii 



= 2Eiiox/NMoj'^ and (y^) = 2Eiioy/NMujy by using the virial relationship ( pT| ) written for the x direction 
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the analogous equations holding for y and z. One then finds 
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where we have explicitly indicated the dependence of E'int , -Ekin x and £^kin y on the trap deformation. By expanding 
the above equations with respect to e one can easily determine the limit entering Eq. dlq) getting, to the lowest order, 



^Qi -1 = v2tJho hm -- 
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with -EkinajC^) — {dEkinx/de)e=o and analogously for y. In order to estimate the diflFerence Ey^^^{Q) — E'^^^^ JQ) 
entering the last equation, we can use the results reported in the next Appendix |D| for the different contributions to 
the kinetic energy in Thomas-Fermi regime. By setting A — LOyju^ and using the relationship e ~ (1 — A) holding for 
small £, one finds that 



^ki„.(0) - i?ki„,(0) = ^i^kin - iV^ 
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which combined with Eq. (C8) finally gives 
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It is worth pointing out that, with logarithmic accuracy, Eq. (CIO) coincides with the result (C4) derived for lu^^i 



APPENDIX D: KINETIC ENERGY IN THOMAS-FERMI REGIME 

By explicitly calculating the various contributions i^kin ± and E'kin z of the total kinetic energy using the method 
developed in Ref. (l7| we find the following results: 
Disk shaped trap, A > 1 
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Cigar shaped trap, A < 1 
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In Eqs. (plD -(D4) 5q represents the surface thickness at z = 0. Its dependence on the relevant parameters of the 
system is [|17| 



1.776/ iVa 



-1/15 
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The above results are consistent with the results of Ref. ||l^ for the total kinetic energy E^kin — -Ekin i + ^kin z , which 
reduces to (17) in the isotropic case A = 1. 
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